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Aberrations of protein-coding genes are a focus of cancer genomics; however, the impact of oncogenes on expression of the 
-50% of transcripts without protein-coding potential, including long noncoding RNAs [IncRNAs], has been largely 
uncharacterized. Activating mutations in the 5Mf oncogene are present in >70% of melanomas, 90% of which produce active 
mutant BRAF^^^ protein. To define the impacts of oncogenic BRAF on the melanocyte transcriptome, massively parallel 
cDNA sequencing [RNA-seq] was performed on genetically matched normal human melanocytes with and without BRAF^^^^^ 
expression. To enhance potential disease relevance by verifying expression of altered genes in 5Mf-driven cancer tissue, 
parallel RNA-seq was also undertaken of two ^Mf^'^^^^-mutant human melanomas. BRAF^^°°^ regulated expression of 1027 
protein-coding transcripts and 39 annotated IncRNAs, as well as 70 unannotated, potentially novel, intergenic transcripts. 
These transcripts display both tissue-specific and multi-tissue expression profiles and harbor distinctive regulatory chromatin 
marks and transcription factor binding sites indicative of active transcription. Coding potential analysis of the 70 unannotated 
transcripts suggested that most may represent newly identified IncRNAs. BRAF-regulated IncRNA 1 [BANCR] was identified as 
a recurrently overexpressed, previously unannotated 693-bp transcript on chromosome 9 with a potential functional role in 
melanoma cell migration. BANCR knockdown reduced melanoma cell migration, and this could be rescued by the chemokine 
CXCLll. Combining RNA-seq of oncogene-expressing normal cells with RNA-seq of their corresponding human cancers may 
represent a useful approach to discover new oncogene-regulated RNA transcripts of potential clinical relevance in cancer. 

[Supplemental material is available for this article.] 



The RAS/RAF/MAPK is pathway is hyperactive in —30% of human 
cancers, and activating mutations in key members of this pathway 
serve as driver mutations in many malignancies (Dankort et al. 
2009; Dhomen et al. 2009; Michailidou et al. 2009). Activating 
mutations in the BRAF oncogene are present in >70% of melanomas, 
-90% of which are BFiAF^^^^ (Davies et al. 2002; Rubinstein et al. 
2010). The recent clinical use of selective BRAF^^^^^ inhibitors 
demonstrates that impairing BRAF^^^^^ signaling evokes tumor re- 
gression in patients with 5iMF-mutant metastatic melanoma (Bollag 
et al. 2010; Flaherty et al. 2010). However, patients eventually relapse 
and succumb to acquired chemoresistance via mechanisms that 
reactivate MAPK signaling (Johannessen et al. 2010; Nazarian et al. 
2010; Villanueva et al. 2010; Wagle et al. 2011). It is therefore clear 
that a greater understanding of oncogenic 5iMF-driven impacts is 
required to yield alternative and complementary therapeutic tar- 
gets. Recent DNA sequencing efforts have focused on protein 
coding exons to provide further insight into genomic mutations in 
melanoma (Berger et al. 2010; Johannessen et al. 2010; Pleasance 
et al. 2010; Prickett et al. 2011; Wei et al. 2011); however, areas of 
the genome that do not encode functional proteins are largely 
unexplored. 

Large-scale transcriptome analyses indicate that -50% of 
transcribed RNAs have no protein-coding potential (Claverie 2005; 
Kapranov et al. 2007; Mercer et al. 2009); however, the functional 
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importance of untranslated RNA in cancer is not fully understood. 
Among these transcripts, long noncoding RNAs (IncRNAs) that 
vary in length from 200 bp to tens of kilobases are a particularly 
important class of noncoding RNA (ncRNA) that have recently 
been shown to regulate diverse functions, including X-chromo- 
some silencing (Borsani et al. 1991; Lee et al. 1999), pluripotency 
(Loewer et al. 2010), epigenetic gene regulation (Rinn et al. 2007; 
Tsai et al. 2010), TP53 (also known as p53) functions (Huarte et al. 
2010), and the DNA damage response (Hung et al. 2011). The 
regulation and function of IncRNAs in cancer is largely unchar- 
acterized, and whether expression of IncRNAs is regulated by 
oncogenic BRAF^^^^^ is unknown. 

We performed RNA-seq on both normal melanocytes (± ec- 
topic BRAF^^^°^) and BRAF-mutdnt human melanomas to define 
the impact of oncogenic BRAF expression on the melanocyte 
transcriptome. BRAF^^°^^ regulated expression of 1027 protein- 
coding transcripts, 39 annotated IncRNAs, and 70 previously un- 
annotated transcripts. These transcripts display both tissue-specific 
and multi-tissue expression profiles and harbor distinctive regu- 
latory chromatin marks indicative of active transcription. We 
validated expression of a subset of these transcripts and also con- 
firmed similar expression patterns in another publicly available 
melanoma RNA-seq data set. One of the novel transcripts most 
highly induced by oncogenic BRAF, which is also recurrently 
overexpressed in melanoma, is a putative IncRNA we named BRAF- 
activated non-coding RNA (BANCR). BANCR regulates a set of 
genes involved in cell migration, including the chemokine 
CXCLll, and BANCR is required for full migratory capacity of 
melanoma cells. CXCLll rescued defective migration in BANCR- 
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knockdown cells, suggesting that CXCLll is an important target 
of BANCR involved in cell migration. These data present a high- 
resolution profile of an oncogene-remodeled transcriptome and 
identify oncogene-regulated IncRNAs and novel RNA transcripts 
of potential interest for future study in human cancer. 

Results 

To assess transcriptome-wide changes induced by BRAF^^^^^, we 
performed RNA-seq on genetically matched primary human me- 
lanocytes after transduction with a lentivirus expressing either 
mutant oncogenic BRAF^^^^^ or red fluorescent protein (RFP) 
marker control. As an additional filter for potential clinical rele- 
vance, RNA-seq was also performed on two primary human mela- 
nomas expressing BRAF^^°°^ (Fig. 1; Supplemental Fig. 1). Transcript 
assembly of mapped reads was performed using two methods 
(Fig. lA): assembly to a reference transcriptome (RefSeq/Gencode) 



using Cufflinks (Trapnell et al. 2010), or unbiased, de novo as- 
sembly using Scripture (Guttman et al. 2010). Assembly to RefSeq 
yielded expression data for annotated protein-coding transcripts 
and RefSeq-annotated ncRNAs. Sequencing reads were also 
mapped to the Gencode reference since it contains additional, 
curated ncRNA annotations that are absent in RefSeq. Addi- 
tionally, all annotated transcripts were deducted from Scripture- 
assembled transcripts to yield a population of novel de novo- 
assembled transcripts (Fig. lA). This approach permitted analysis 
of BRAF-dependent regulation of protein-coding transcripts, 
annotated IncRNAs, and novel transcripts. BRAF^^°^^ concor- 
dantly altered expression twofold or more of 1027 protein-cod- 
ing transcripts, 39 annotated IncRNAs, and 70 novel transcripts 
(Fig. IB). Unsupervised hierarchical clustering in both melano- 
cytes expressing BRAF^^°°^ as well as 5iMF^^^^^-positive human 
melanoma tissue demonstrated concordant expression changes 
for the majority of protein coding transcripts (Fig. IC; Supple- 
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Figure 1. Oncogenic BRAF^^°°^ remodels the transcriptome in melanoma. (A) Schematic of experimental workflow and RNA-seq data analysis. (6) Pie 
chart showing number and categories of transcripts regulated by BRAF^^°°^ Numbers of transcripts changed in all three samples in the same direction ± 
twofold versus pooled normal melanocyte control. (C) Heatmap showing hierarchical clustering of protein-coding transcript expression in primary 
melanocytes overexpressing BRAF^^°°^ and in primary melanomas (± twofold, changed in all three samples). (D) Top 10 unique gene ontology (GO) 
terms associated with protein-coding transcripts up-regulated by BRAF^^°°^. 
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mental Table 1). The top 10 unique gene ontology (GO) terms 
associated with BRAF up-regulated protein-coding transcripts 
agree with GO terms previously associated with 5i^^F-mutant 
melanoma (Fig. ID; Kannengiesser et al. 2008; Packer et al. 2009), 
supporting the potential clinical relevance of the RNA-seq data 
obtained. 

In normal human melanocytes, BRAF^^^^^ significantly al- 
tered expression of 39 annotated IncRNAs and 70 novel tran- 
scripts that also exhibited twofold or more altered expression in 
BRAF^^^^^-expressing primary melanomas (Fig. IB; Supplemental 
Tables 2, 3). We validated regulation by BRAF for 12/12 annotated 
transcripts tested and 13/13 novel transcripts tested (Supplemental 
Tables 4, 5, respectively). We then assessed expression levels of 
these 109 transcripts in multiple tissues using publicly available 
ENCODE RNA-seq data and found that transcripts displayed both 
tissue-specific and multi-tissue expression profiles (Fig. 2A; Sup- 
plemental Table 6). DNase 1 hypersensitive sites (DNaseH) are 
indicative of open chromatin and are commonly enriched at 
transcription start sites (TSSs) and enhancers. We interrogated 
ENCODE data for DNaseH sites in primary human melanocytes 
and found enrichment at the TSSs of these transcripts (Fig. 2B). 
In addition, we also looked for chromatin modifications in other 
ENCODE cell types where these transcripts are most abundant 
since no other ENCODE data exist for primary melanocytes. TSSs 
were enriched for H3K4me3, RNA pol II binding, and DNaseH in 
HIES cells (Fig. 2C-E) and in K562 leukemia cells (Fig. 2F-H). 
Loading of RNA pol II and the presence of chromatin modifi- 
cations associated with active transcription further validate ex- 
istence and expression of these annotated IncRNAs and novel 
transcripts in melanocytes and other tissues. In addition, we also 
analyzed ChlP-seq data for 77 transcription factors (ENCODE) 
and found that most transcripts (79/109) bound at least one 
transcription factor in its promoter region, and positive binding 
was more highly represented in transcripts with broad tissue 
expression (Supplemental Fig. 2; Supplemental Table 7). A rep- 
resentative example of an annotated IncRNA regulated by BRAF 
in melanoma, expressed in other tissues that harbors transcrip- 
tion-associated chromatin marks is shown in Supplemental 
Figure 3. 

Next, we focused on exploring putative transcripts discovered 
by unbiased de novo assembly of RNA-seq data using previously 
validated techniques (Guttman et al. 2010; Trapnell et al. 2010). By 
use of this approach, 70 unannotated transcripts were identified as 
concordantly changed twofold or more in both melanocytes 
overexpressing BRAF^^^^^ and in primary melanoma tissue sam- 
ples (Fig. IB; Supplemental Table 3). To assess the protein-coding 
potential of these transcripts, we used a coding potential calculator 
(CPC) algorithm previously shown to discriminate coding from 
noncoding transcripts with high accuracy (Kong et al. 2007). Posi- 
tive CPC scores indicate positive protein-coding potential, whereas 
negative CPC scores indicate low coding-potential. In agreement 
with this, of 1500 randomly chosen annotated ncRNA transcripts 
analyzed, most yielded negative CPC scores, whereas 1500 ran- 
domly selected protein-coding genes yielded positive CPC scores 
(Fig. 3 A). Newly discovered BRAF-regulated transcripts returned 
a similar CPC profile to the control ncRNA population, suggesting 
that many of these transcripts may be IncRNAs (Fig. 3 A). Next, to 
search for biologically important, highly expressed, BRAF-regu- 
lated IncRNAs, we assessed expression levels of up-regulated novel 
transcripts with negative CPC scores, in a publicly available mel- 
anoma RNA-seq data set published recently (Berger et al. 2010). 
This study did not search for IncRNAs but instead focused on 



identifying novel gene fusions and other events driven by un- 
derlying aberrations in genomic DNA. Re-mapping the RNA-seq 
data from Berger et al. (2010; including intergenic regions) allowed 
us to derive expression values (using Cuffdiff) for the novel IncRNAs 
identified in our study Clustering of these data revealed recurrent, 
overexpression of a novel IncRNA that we named BANCR (Fig. 3B- 
H; Supplemental Fig. 4). 

RNA-seq data histograms and de novo assembly using Scrip- 
ture indicated that BANCR is a four-exon transcript of 693 bp that 
is highly induced by BRAF^^^^^ in melanocytes and is overex- 
pressed in primary tumors (Fig. 3C). We also detected the BANCR 
transcript using the Trinity algorithm (Grabherr et al. 201 1), which 
does not rely on mapping of RNA-seq reads to a reference genome 
sequence. This complementary approach provided orthogonal 
validation of BANCR transcript architecture, and in further support 
of this, there are spliced ESTs that overlap with our BANCR gene 
prediction (e.g., EST W58600). De novo assembled transcripts also 
align faithfully with RNA-seq reads in the independent melanoma 
data set of Berger et al. (2010; Supplemental Fig. 4B). The CPC score 
for BANCR is similar to that derived for the well-characterized 
IncRNAs HOTAIR (Rinn et al. 2007) and XIST (Borsani et al. 1991), 
whereas control protein-coding transcripts such as BRAF and TP53 
yield positive CPC scores (Fig. 3D). BANCR expression levels were 
increased in 4/5 additional BiMF^^^*^^ melanoma samples, to levels 
as high as 20- to 40-fold (Fig. 3E). We also validated up-regulation 
of BANCR by both BRAF^^^°^ and mutant active NRAS in mela- 
nocytes by qPCR (Fig. 3F). Northern blot analysis revealed up- 
regulation of a single band corresponding to BANCR RNA in me- 
lanocytes expressing BRAF^^^^^ and in the BRAF-mutdnt mela- 
noma cell line sk-mel-5 (Fig. 3G). BANCR is transcribed at a locus 
on chromosome 9 and lies in a gene desert 40 kb downstream from 
TJP2 and 20 kb upstream of FAM189A2 (Fig. 3H). Expression of the 
most proximal genes upstream of and downstream from BANCR 
was not significantly changed (Supplemental Table 8), suggesting 
that BANCR is not a c/s-regulator of neighboring genes as has been 
described for a subset of IncRNAs (Pasmant et al. 2007; Yap et al. 
2010). By exploring multi-tissue RNA-seq data from ENCODE, we 
found that BANCR expression is largely melanocyte/melanoma- 
specific (Fig. 2A, marked with an asterisk), but we did detect low 
level expression in K562 leukemia cells and HIES cells (Fig. 2A; 
Supplemental Fig. 5). In addition, the BANCR locus is enriched for 
chromatin marks associated with active transcription, and in- 
terestingly, the BANCR promoter was shown to be enriched for 
STATl transcription factor binding in HeLa-s3 cells (ENCODE 
ChlP-seq data) (Supplemental Fig. 5). 

To explore if BANCR may be relevant to BRAF^^^^^-driven 
oncogenic impacts, we undertook shRNA-mediated BANCR de- 
pletion in melanoma cells. First, gene expression profiling was 
performed using complementary DNA (cDNA) microarrays on 
Colo829 BRAF^^^^^ melanoma cells in which BANCR expression 
had been reduced to <25% using shRNAs (Fig. 4A). Eighty-eight 
genes changed significantly following BANCR knockdown using 
independent BANCR shRNA sequences (Fig. 4B). Genes repressed 
by BANCR depletion were enriched for GO terms, including loco- 
motory behavior and chemotaxis (Fig. 4C), suggesting a potential 
role for BANCR in melanoma cell motility. Consistent with this, 
BANCR depletion impaired migration of melanoma cells without 
affecting viability and proliferation (Fig. 4D,E; Supplemental Fig. 
6A), indicating that BANCR is at least partially necessary for mel- 
anoma cell migration in vitro. We next tested if factors down- 
regulated by BANCR loss were able to rescue defective cell migra- 
tion. We found that the chemokine CXCLll, which is down- 
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Figure 2. Multi-tissue expression profile of annotated IncRNAs and novel transcripts. (A) Publicly available RNA-seq data (ENCODE project) for multiple 
cell types were downloaded from the UCSC Genome Browser. The melanoma column expresses an aggregate from our RNA-seq data indicating if each 
transcript was concordantly up-regulated (yellow) or down-regulated (blue) by BRAF in melanoma. Red/white scaling indicates expression level in other 
cell types. (B) Melanocyte DNaseH (ENCODE) distribution over the 109 BRAF-regulated transcripts is shown. The following charts (C-H) depict signal 
distribution (ENCODE) over the 109 BRAF-regulated transcripts. (C) Distribution of H3K4me3 in HI ES cells. (D) Distribution of RNA pol II in HI ES cells, 
(f ) Distribution of DNaseH in HI ES cells, (f) Distribution of H3K4me3 in K562 cells. (C) Distribution of RNA pol II in K562 cells. (H) Distribution of DNaseH 
in K562 cells. (TSS) Transcription start site; (TES) transcription end site. 
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Figure 3. Oncogenic BRAF^^°°^ regulates expression of many novel transcripts. (A) Frequency distribution plot showing coding potential analysis (CPC 
scores) of novel transcripts discovered by de novo assembly. One thousand five hundred random annotated ncRNAs and 1500 random protein-coding 
genes were also analyzed for reference. (S) Analysis of BRAF up-regulated novel IncRNAs in a publicly available melanoma RNA-seq data set. We remapped 
and reanalyzed Berger et al. (201 0) data and derived expression values for novel transcripts up-regulated in all three of our samples. Robustly up-regulated 
novel transcripts (with negative CPC scores) were clustered with Berger et al. (201 0) data revealing recurrently highly expressed IncRNAs. (C) Histograms 
of raw RNA-seq data in control sample, in melanocytes overexpressing BRAF^^°°^ and in both primary melanomas (Mel. 1, Mel. 2). Scripture and Trinity 
assembly for BANCR is also shown, /-axis is number of RNA-seq reads normalized for mapping variation. (D) CPC score for BANCR. Examples of IncRNAs 
HOTAIR and X/57and protein-coding genes BRAF and TPS 3 are shown for reference, (f ) BANCR expression in primary melanomas validated by qRT-PCR. (f ) 
6/4NC/? expression in melanocytes overexpressing BRAF^^°°^ or NRAS*"^^^ measured by qRT-PCR. Data are means from three experiments ± SD. (C) BANCR 
Northern blot using RNA from control melanocytes (MCs), melanocytes overexpressing BFIAF^^^^^, or sk-mel-5 melanoma cells. (H) Schematic of BANCR locus. 



regulated by 5ANCi^ knockdown (Supplemental Fig. 6B), was able 
to rescue the migration defect without significantly altering 
baseline migration (Fig. 4F,G; Supplemental Fig. 6C). These data 
indicate that BANCR is a novel BRAF^^^^^-induced IncRNA that 
regulates melanoma cell migration in vitro in part by regulating 
expression of CXCLll. 



Discussion 

In this study we performed RNA-seq of BRAF^^^^^-expressing 
primary human melanocytes to characterize how this oncogene 
remodels the melanocyte transcriptome. To enhance cancer rele- 
vance, these data were filtered through RNA-seq data from duplicate 
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Figure 4. BANCR loss impairs melanoma cell migration. (A) Quantification of BANCR knockdown by qPCR using two independent shRNAs. (6) 
Microarray analysis of Colo829 melanoma cells following shRNA knockdown of BANCR. (C) Significant gene ontology (GO) terms associated with genes 
repressed by BANCR loss. (D) Representative images showing amount of cell migration at day 0 and day 2. Melanoma cells were infected with lentivirus 
expressing either nontargeting control shRNA (control) or one of two duplicate independent shRNA sequences targeting BANCR (shRNA 1, shRNA 2). 
Green highlight, shown for clarity, is a mask applied by the analysis software for accurate, unbiased migration quantification, (f ) Quantification of replicate 
images (n = 6 per timepoint) taken over migration assay timecourse. Values are means from two independent experiments ± SD. Statistical analysis performed 
by one-way AN OVA, **P< 0.01, *P< 0.05 versus control, (f) Quantification of replicate images (n= 3 per timepoint) taken over migration assay timecourse ± 
CXCL1 1 . Cells were infected with shRNA control lentivirus and treated with vehicle (PBS), infected with BANCR shRNA, and treated with vehicle or infected with 
BANCR shRNA and treated with CXCLl 1 (10 ng/mL). (C) Representative images showing amount of cell migration at day 0 and day 2. 



BRAF^^^^^ mutant human melanomas to characterize BRAF^^°°^- 
regulated RNA transcripts that are also recurrently altered in mel- 
anoma. By use of this approach, we identified 39 annotated 
IncRNAs regulated by BRAF^^°°^ that are also concordantly 
expressed in melanoma. By using de novo transcript assembly 
we identified 70 additional un-annotated BRAF^^°^^-regulated 
transcripts that were also concordantly expressed in melanoma. By 
comparing these data to a reanalyzed, publicly available mela- 
noma RNA-seq data set, we identified a list of priority BRAF target 
transcripts that are recurrently altered in cancer. From these tran- 
scripts, the putative IncRNA, BANCR, was identified as a potential 
regulator of melanoma cell migration. 



Large-scale transcriptome analyses indicate that —50% of 
transcribed RNA has no protein-coding potential (Claverie 2005; 
Kapranov et al. 2007; Mercer et al. 2009). How this "dark matter" of 
the genome is altered in cancer, and the functional importance of 
such alterations, is only beginning to be explored. W^hereas ncRNA 
species such as microRNAs are established as being important in 
cancer, more recently described IncRNAs are less well studied 
in this context. Altered expression of IncRNAs has been reported 
in several cancers: For example, HOTAIR is up-regulated in primary 
and metastatic breast tumors, and high levels of HOTAIR were as- 
sociated with metastasis and poor survival (Gupta et al. 2010); 
lincRNAp21 is regulated by TP53 and mediates the TP53 response 
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(Huarte et al. 2010); SPRY4-IT1 is overexpressed in melanoma 
(Khaitan et al. 201 1); the IncRNA PANDA is aberrantly expressed in 
breast cancer and is implicated in regulating the DNA damage re- 
sponse (Hung et al. 2011); MALAT-1 is associated with metastasis 
in non-small-cell lung cancer (Ji et al. 2003); and ANRIL is involved 
in epigenetic silencing of pl5 in leukemia (Yu et al. 2008). 
Whether altered IncRNA expression represents noise resulting 
from tumor heterogeneity and genomic instability or whether 
this is driven by oncogenic driver events has not been determined 
systematically. Employing RNA-seq and using BRAF^^^^^ mela- 
noma as a system to interrogate oncogene-regulated IncRNAs 
allowed visualization as to how an oncogene transforms the 
transcriptional landscape in normal cells. It also facilitated the 
identification of several novel transcripts and IncRNAs that are 
frequently dysregulated in melanoma, many of which are also 
expressed in other tissues. 

Recent RNA sequencing of eight short-term melanoma 
cultures originating from stage 4 tumors and two melanoma cell 
lines identified 11 novel melanoma gene fusions and 12 novel 
read-through transcripts (Berger et al. 2010). This insightful 
study, however, did not seek to discover novel RNA transcripts, 
which may have been technically challenging due to the con- 
siderably lower sequencing depth compared with that of 
our present study. Berger et al. (2010) obtained an average 
of 14.2 million total reads and 8.4 million mapped reads per 
sample, whereas our present study achieved an average of 151 
million total reads and 132 million mapped reads per sample 
representing 10X-15X increased coverage (Supplemental Table 
9). We subsequently remapped and reassembled Berger et al. 
(2010) data and included intergenic regions and our de novo 
assembled novel transcripts, which increased mapped reads 
from 8.4 million to around 20 million (Supplemental Table 10). 
This allowed interrogation of the expression of the newly 
discovered novel transcripts described here and revealed that 
BANCR is recurrently overexpressed in melanoma. Since BANCR 
is overexpressed in both our primary melanomas and metastatic 
samples from Berger et al. (2010), BANCR may play a potential 
role throughout the spectrum of disease progression. Consistent 
with this, knockdown of BANCR in metastatic melanoma cell 
lines inhibited melanoma cell migration (Fig. 4). 

Our findings that CXCLll was down-regulated by BANCR 
loss in migration-deficient melanoma cells and that migration 
could be rescued by reintroduction of recombinant CXCLll are 
consistent with previous findings. For example, it has been shown 
that the anti-tumor drugs dimethylfumarate (DMF) and dacarba- 
zine (DTIC) reduced migration and CXCLll expression in mela- 
noma cells in vitro (Valero et al. 2010). In addition, melanoma cells 
lacking the CXCLll receptor CXCR3 exhibited reduced migration 
in a murine model of melanoma metastasis (Kawada et al. 2004). 
Our data further implicate CXCLll as being an important medi- 
ator of melanoma cell migration and identify BANCR as an im- 
portant regulator. 

In conclusion, we have used RNA-seq to identify annotated 
and novel transcripts regulated by oncogenic BRAF^^°°^. In- 
terrogation of novel IncRNAs up-regulated by BRAF^^°°^ that were 
also recurrently up-regulated in melanoma identified BANCR, a 
novel IncRNA that regulates melanoma cell migration. The com- 
bination of RNA-seq of normal cells expressing relevant oncogenes 
with RNA-seq of the corresponding oncogene-expressing human 
cancers may provide a useful general approach to discover new 
oncogene-regulated RNA transcripts of possible clinical relevance 
to human cancer. 



Methods 

Cell culture 

Primary human melanocytes were isolated from fresh, discarded 
surgical foreskin specimens obtained from Lucile Packard Chil- 
dren's Hospital at Stanford through and Institutional Review 
Board-approved protocol. Isolation performed as previously de- 
scribed (Chudnovsky et al. 2005). Primary melanocytes were cul- 
tured in melanocyte medium 254 (Invitrogen). 293Tand sk-mel-5 
cells were cultured in DMEM + 10% FCS (Invitrogen). All cells were 
cultured in a 37°C, 5% CO2 humidified incubator. Colo-829 cells 
were cultured as above using RPMI + 10% FCS (Invitrogen). 

Lentiviral infections and shRNA vector construction 

gf^pvsooE NRAS'^^^^ were cloned into the lentiviral vector 
pRRL-sin-cPPT-hPGK-WPRE. Helper packaging plasmids used were 
pucMDG, and pCMVA8.91. 293T cells were transfected overnight 
with virus production vectors using calcium phosphate trans- 
fection, and the cell medium was replaced with DMEM +10% FCS 
the following day. Transfected cells were transferred to a 32°C, 5% 
CO2 incubator overnight for virus production. Melanocytes were 
infected with diluted high titer lentivirus by centrifugation in 
10-cm culture plates at 32°C in the presence of 5 jxg/mL polybrene 
(Sigma). As a control, a population of normal melanocytes pooled 
from three donors was infected with virus transmitting RFP. Cus- 
tom-made shRNAs were cloned into the lentiviral vector pGIPZ 
(Openbiosystems) and were designed using RNAi Central shRNA 
des. tool (http://cancan.cshl.edu/RNAi_central/RNAi.cgi?type= 
shRNA). Cells were infected twice, 24 h apart as above to achieve 
best knockdown. All shRNA sequences are shown in Supple- 
mental Table 2. 

Western blotting 

Cells were lysed in cell lysis buffer (no. 9803, Cell Signaling); 
protein concentration was determined using Bio-Rad protein assay 
dye (no. 500-0006, BioRad); and polyacrylaminde gel electropho- 
resis was performed using NuPAGE 10% Bis-Tris gels (Invitrogen). 
Proteins were transferred to a Amersham Hybond-P PVDF mem- 
brane (GE Healthcare). Antibodies used were as follows: BRAF (sc- 
5284, Santa Cruz Biotechnology), pERK (no. 4377, Cell Signaling), 
total ERK (no. 9102, Cell Signaling), beta-actin (A2668, Sigma), 
anti-Mouse HRP (NA931), and anti-Rabbit HRP (NA934) secondary 
antibodies (GE healthcare). Membranes were developed using ECL 
reagent (GE Healthcare). 

Northern blotting 

RNA was separated by RNA gel electrophoreses and transferred to 
Hybond-N+ nylon membrane (Amersham) overnight using 20 X 
SSC buffer. The membrane was cross-linked using a UV Stra- 
talinker. The Northern probe was made by PCR-cloning BANCR 
cDNA followed by labeling with a-^^P-labeled dCTP (Perkin 
Elmer). Northern membrane and probe were incubated overnight 
at 60°C in a rotating incubator. The membrane was then washed 
with 2X SSC buffer +0.1% SDS and exposed overnight at -80°C for 
visualization. 

RNA extraction, cDNA library construction, 
and cDNA sequencing 

Total RNA was isolated with an RNeasy mini kit according to the 
manufacturer's instructions (Qiagen). cDNA sequencing libraries 
were prepared using the Illumina mRNA sequencing prep kit (RS- 
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100-0801), which selects for poly-A containing mRNAs. Two mi- 
crograms of total RNA from each sample was used for library prep- 
aration. High-throughput, 101 -bp paired-end cDNA sequencing 
was performed by Elimbio. 

RNA-seq data analysis 

Sequencing reads were mapped to NCBI build 36.1/hgl8 using 
"TopHat," and mapped sequences were aligned to RefSeq and 
Gencode transcriptomes using "Cufflinks" (Trapnell et al. 2010). 
Transcripts were also assembled using "Scripture," which performs 
unbiased de novo transcript assembly without using a reference 
transcriptome (Guttman et al. 2010). For fold change calculations, 
RNA-seq data were compared with data from a population of nor- 
mal melanocytes pooled from three donors, infected with Anti- 
virus transmitting RFR Transcript coding potential calculation was 
performed using "CPC calculator" (Kong et al. 2007). GO terms 
associated with significantly changed protein-coding genes were 
calculated using DAVID (Huang da et al. 2009). 

ENCODE data analysis 

For multi-tissue expression analysis, paired-end RNA-seq data for 
HIESC, K562, HSMM, NHLF, HUVEC, GM12878, HELA, HEPG2, 
MCF7, and NHEK cells (all ENCODE project) were downloaded 
from the UCSC Genome Browser. Expression was quantified with 
Cufflinks using a reference transcriptome composed of RefSeq, 
GENCODE, and predicted Scripture-derived transcripts. FPKM 
values were hierarchically clustered and visualized using Java 
Treeview. 

Transcription factor binding 

Binding sites for 77 transcription factors were obtained from the 
UCSC table browser (ENCODE transcription factor ChlP-seq peak 
files). Transcript promoter regions were defined as 3 kb upstream of 
to 1 kb downstream from the TSS. With this, we generated a regu- 
latory matrix of transcripts as rows and 77 transcription factors as 
columns. A 0 indicates no binding; 1, positive binding detected. 

ChlP-seq data for K562 and HIES cells from the ENCODE 
project (Birney et al. 2007) for H3K4me3 and RNA pol II (signifi- 
cant called ChlP-seq peaks) were downloaded from the UCSC 
Genome Browser. We also downloaded data for DNasel hyper- 
sensitive regions for primary human melanocytes, K562, and HIES 
cells. Wiggle files were converted to bedGraph format, and the 
average signal was calculated beginning upstream of the TSS by 
a distance of 50% of the gene body (TSS to transcription stop site) 
and continuing downstream from the transcript by the same 50% 
distance. These data were sorted into 100 bins normalized to the 
transcript length. 

Reverse transcription real-time qPCR 

RNA was reverse transcribed using the iScript Reverse transcription 
kit (BioRad). qPCR was performed using the Maxima SYBR Green/ 
ROX qPCR master mix (Thermo Scientific). Primers for qPCR were 
designed using "OligoPerfect" primer design software (Invitrogen), 
and the primer specificity and efficiency/dynamic range was con- 
firmed to be linear over four orders of magnitude during optimi- 
zation prior to use in experiments. All qPCR primers used are listed 
in Supplemental Table 1. 

Microarray analysis 

Microarrays were performed using GeneChip Human Genome 
U133 Plus 2.0 array chips (Affymetrix) on RNA extracted from 



Colo829 melanoma cells ± BANCR knockdown. cDNA synthesis, 
labeling, and array procedures were conducted at the Stanford PAN 
facility. 

Melanoma RNA 

RNA from primary human melanomas was purchased from 
Asterand. H and E histological images of melanoma samples and 
tumor grading information was also provided by Asterand. 

BRAF and NRAS genotyping 

The BRAF^^^^ locus and NRAS^^^ locus were amplified and purified 
from melanoma cDNA. BRAF locus primers were as follows: for- 
ward, agcttgtatcaccatctccatatca; reverse, accataaaaattatctggtccctgt. 
NRAS locus primers were as follows: forward, ctgagtacaaactggtg 
gtggtt; reverse, ggcttgttttgtatcaactgtcc. Products (BRAF 401 bp, 
NRAS 386 bp) were resolved by gel electrophoresis, purified, and 
sequenced by Sanger sequencing. 

Migration assays 

A modified "scratch" migration assay that does not involve vari- 
able, physical wounding of cell cultures was used to assess cell mi- 
gration (Ibidi, see http://www.ibidi.de). Cells were seeded at 40,000 
cells per chamber and allowed to attach overnight. The following 
day, culture inserts were removed and light microscopy images ac- 
quired (three per condition). Cells were cultured during the migra- 
tion assay under normal conditions, and images were subsequently 
acquired 24 h and 48 h later. Images were analyzed using automated 
image analysis software (http://ibidi.wimasis.com/). For CXCLll 
experiments, CXCLll (R&D Systems) diluted in PBS was applied at 
10 ng/mL. 

Data access 

All RNA-seq data and cDNA microarray data have been submitted 
to the NCBI Gene Expression Omnibus (GEO) (http://www.ncbi. 
nlm.nih.gov/geo/) under accession number GSE33092. 
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